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Abstract. In this study, after we have briefly introduced the standard Gross- 
Pitaevskii equation, we have suggested fractional Gross-Pitaevskii equations to 
j^ \ investigate the time-dependent ground state dynamics of the Bose-Einstein 

^ ' condensation of weakly interacting bosonic particle system which can includes non- 

Markovian processes or non-Gaussian distributions and long-range interactions. Only 
we focused the time-fractional Gross-Pitaevskii equation and have obtained solutions 
1,^ of the standard Gross-Pitaevskii and time-fractional Gross-Pitaevskii equations for 

*^ , attractive and repulsive interactions in the case external trap potentials V{x) = and 

lO ' optical lattice potential V {x) — ± sin'^ x by using Homotopy Perturbation Method. We 

^^ ■ have found that the Homotopy Perturbation Method solutions of the Gross-Pitaevskii 

equation for these potentials and interactions are the same analytical results of it. 
Furthermore we have also found that solutions of the time-fractional Gross-Pitaevskii 



m 

fvq , equation for these potentials and interactions can be given in terms of Mittag-Leffler 

function. The solutions of the time-fractional Gross-Pitaevskii equation provide that 
the time evolution of the ground state dynamics of Bose-Einstein condensation of 
bosonic particles deviates exponential form, and evolutes with time as stretched 
exponentially. 



PACS numbers: 67.85.Hj, 67.85.Jk, 02.30.Jr 
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1. Introduction 

The zero temperature Bose-Einstein condensation (BEC) of dilute or weakly interacting 
bosons are generally described by Gross-Pitaevskii equation which is s a self-consistent 
mean field nonlinear Schrodinger equation [HE]- It is clearly known that this 
equation does not include memory effects of the non-Markovian processes, non-Gaussian 
distribution of particles, the long-range inter-particle interaction effects or the fractal 
structure of the interacting space between particles in condensate phase. However, it is 
shown that some experimental observations in a real real dilute or weakly interacting 
bosonic systems deviates from theoretical predictions i.e. some experimental results 
do not comply with theoretical calculations. For example, for dilute atomic gases 
that involve interacting atoms or molecules such as ^^Rb, ^^Na and ^Li which are 
trapped in a harmonic oscillator potential, there are significant differences between the 
theoretical and experimental transition temperatures in Bose-Einstein condensation. 
Similarly, for ^He, the experimentally observed value of the transition temperature to 
BEC is Tc = 2.17 K, whereas the calculated one is T^ = 3.10 K in conventional Bose- 
Einstein thermostatistics [H |5l El El [HI |9]. These contractions between experimental 
observation and theoretical predictions shows that the memory effects of the non- 
Markovian processes, non-Gaussian distribution of particles, the long-range inter- 
particle interaction effects or the fractal structure of the interacting space in real bosonic 
systems have a predominant role below the critical temperature Tc [3]. For this reason, 
we can claim that this standard approaches remain insufficient in the investigation 
ground state dynamics of the physical behavior of weakly interacting quantum boson 
gases at low temperatures i.e. below the T^. At finite temperature below the Tc and even 
though at zero temperature, the non-Markovian processes, non-Gaussian distribution 
of particles, the long-range effects between particles or the fractal structure of the 
interacting space in the Bose-Einstein condensate phase must be taken into account 
in the model. However, to the best of our knowledge on the literature, there does not 
exist a sound mathematical basis for these efforts in the literature. 

In this study, we will adopt fractional calculus, which involves fractional derivatives 
and integrals. With the help of fractional calculus, it is possible to take into account 
memory effects of the non-Markovian processes, non-Gaussian distribution of particles, 
the long-range inter-particle interaction effects or the fractal structure of the interacting 
space [3l[ini[II]. Because in the fractional calculus concept, these unusual effects that are 
disappearing in the Markovian process can be represented the space-fractional and time- 
fractional derivative operators [12l[13]. Therefore, here we will define fractional Gross- 
Pitaevskii equations, and we will focus only time-fractional Gross-Pitaevskii equation to 
investigate ground state dynamics of the Bose-Einstein condensate of dilute or weakly 
interacting boson gas has repulsive or attractive interactions for zero external potential 
and optical lattice potentials at zero temperature. 

As it is known that time- dependent Gross-Pitaevskii equation is a nonlinear 
partial differential equation and it is difficult to find analytical solution of this type of 
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equation for especially complex potentials. Therefore, the numerical solution methods 
are generally used to find solution of them. In order to solve fractional ordinary 
differential equation, integral equation and fractional partial differential equations in 
physics and other areas of science, several analytical and numerical methods have 
been proposed. The most commonly used ones are; Adomain Decomposition Method 
(ADM) [H Ha Hg [m [m Hg, variational Iteration Method (VIM) [IH HH [191 [201 |2l] , 
Fractional Difference Method (FDM) |22], Differential Transform Method (DTM) [23] , 
Homotopy perturbation Method (HPM) pll [25] . Also there are some classical solution 
techniques e.g. Laplace transform method. Fractional Green's function method, Mellin 
transform method and method of orthogonal polynomials [22]. In this study, we 
will use the HPM method to solve fractional Gross-Pitaevskii equations. HPM has 
been proposed by He [26l [27| to solve linear and nonlinear differential and integral 
equations. This method, which is a coupling of the traditional perturbation method 
and homotopy in topology, deform continuously to simple problem which easily solved. 
The HPM can be easily applied to Volterra's integro-differential equation [28], to 
nonlinear oscillators [29], bifurcation of nonlinear problems [30], bifurcation of delay- 
differential equations [31], nonlinear wave equations [32], boundary value problems 
[33] . quadratic Ricatti differential equation of fractional order [2l] and to other fields 
[3l[35l[36l[371[38l[39l[lQl[lIl[l2]. This HPM yield very convergence of the solution 
series in most cases, usually only a few iteration leading to very accurate solutions. 

This paper is organized as follows. In Section 2 we briefly review Gross-Pitaevskii 
equation. In Section 3, we present the stationary solution of the Gross-Pitaevskii 
equation. In Section 4, we define Gross-Pitaevskii equation with fractional orders. In 
Section 5, we summarize basic definition of the fractional calculus. In Section 6, the 
Homotopy perturbation method is introduced. In Section 7, we discuss several example. 
Finally Section 8 is devoted to conclusions. 

2. Gross-Pitaevskii equation 

The properties of a Bose-Einstein condensate at zero temperature are well described 
by the macroscopic wave-function ip {x, t) whose evolution is governed by the Gross- 
Pitaevskii equation [H [2] which is a self-consistent mean field nonlinear Schrodinger 
equation: 

ih^i,{x,t) = u|-v' + \/(x)+^|^(x,t)'|U(x,t) , xen<zR\i) 

V' (x, t) = , xeV = dVL, (2) 

if) {x, 0) = tpQ{x) , X e fi; (3) 

where g, m and V (x) are the interacting parameter between particles, mass of the 
particles and external potential applying to the particle systems, respectively. The 
interacting parameter i.e. coupling constant g is defined as g = Airhas/m where a^ 
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the scattering length of two interacting bosons. Couphng constant determines the 
interaction types between particles. For g < and g > interactions between bosonic 
particles are attractive and repulsive, respectively. 

It is known that there are two important invariant of Eq. (1) which are the 
normalization of the wave-function 

N{iIj)= [ \'^{x,t)\^dx = N{iPo)= [ \'^o{x)\^dx = l , t>0 (4) 

Jn Jn 

and the energy functional for t > 

^ W = / [^ IVV^ ix,t)f + Vix) \^ix,t)f + I \^Pix,t)\^] dx ^E{i;o) (S; 
Jn iz / . 

where A^ (ip) is particle number and E {ip) is the energy of the particle systems in 

condensate phase. It can be seen from Eq. (5) the energy functional E (ip) consist of the 

three parts i.e., kinetic energy, potential energy and interacting energy. 

The magnitude square of the eigenfunction, {ip {x,t)\ , represents the probability 

density of finding a particle at position x and time t. Stationary state of the Bose- 

Einstein condensed system is the independent of the time. To find stationary solution 

of Eq. (1), we write the wave function as 

ij {x, t)=i) (x) e-'^''^ (6) 

where ^ is the chemical potential of the condensate and ip {x) is a function independent 
of time. Substituting (6) into (1) yield the equation 



) = (^--V^ + V{x)+g\ip {xf'\\ V^ (x) , xeVL(ZR'^ 



Hip (x) = {--y~ + V [X) + g 
for ip [x) under the normalization condition 

\ip (x)| dx = 1 



(7) 



(8) 

On the other hand, any eigenvalue /i can be computed from its corresponding 
eigenfunction ip by 



fi = fi{tp) 



I \ViP {x)f + V (x) \ip {x)f + g\ij (x)l' 



dx . 



(9) 



3. Soliton or stationary solutions of the Gross-Pitaevskii equation 

The most simple solution of the Eq. (1) is the stationary solution of its for V" (x) = 0. 
When d = 1, V (x) = the soliton solutions can be obtained (See j2]). For repulsive 
interactions i.e. g > the solution of Eq. (1) is given by 



ipd (x) = tpd {x, 0) 



■ tan h (y/JIx) 



(10) 



This solution is known as dark solution. Hence we can write the time-dependent dark 
soliton as 



ipd{x,t) 



tan/i {y/Jlx) e 



-ifit/h 



:ir 
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which corresponds to a local minimum of the density distribution \ilj{x,t)\ . On the 
other hand, for attractive interactions i.e. g < 0, the solution of the Eq. (1) is given by 

ipb{x) =il)b{x, 0) = J -sec h{y^x) . (12) 

This solution is known as bright solution. Similarly, time-dependent bright soliton 
solution is given by 

i^b {x, t) = ,p sec h i^x) e-'^'l^ (13) 

which corresponding to a local maximum of the density distribution \il} (x, t)| . These 
solutions can be generalized easily in the progressive wave form. Generalized dark soliton 
solution of Eq. (1) can be written in the one-dimensional form 

^ (x, t)^ = J- tan h[^{x-vt- xq)] e'^^*/^ (14) 

and similarly generalized bright solution solution of Eq. (1) is given by 

ip (x, t)^ = J- sec h[^{x-vt- xo)] e''^''^ (15) 

where v is the speed of the soliton. Also a;o and to are respectively initial position and 
time. Solitons are quantized vertex in Bose-Einstein condensate system, which can be 
observation experimentally. The shape of the solitons are change depending on external 
potential. We here discussed only V" (x) = case for simplicity. 

4. Gross-Pitaevskii equations with fractional derivatives 

We remark in introduction that the memory effects of the non-Markovian processes, non- 
Gaussian distribution of particles, the long-range inter-particle interaction effects or the 
fractal structure of the interacting in a bosonic system below the condensate temperature 
Tc space can leads to fractional dynamics. Simply we say that all non-Markovian 
stochastic processes of the particles or waves in condensate phase of bosonic system 
generate the time-fractional Gross-Pitaevskii equation ignoring boundary conditions and 
sources. Similarly, non-Gaussian distribution of the particles or waves in condensate 
phase leads to space-fractional Gross-Pitaevskii equation. On the other hand the 
combination of the non-Markovian and non-Gaussian behavior yield time and space- 
fractional Gross-Pitaevskii equation. These equation are defined in Ref. [3]. In the light 
of these knowledge, the one-dimensional time-fractional Gross-Pitaevskii equation can 
be written in the form 

ih^ (x, t) = ,D]-^ i-^^'' + V{x)+g\ilj (x, tf\\ ^ (x, t) (16) 

where oDl~°' is the Riemann-Liouville fractional integral operator. Similarly the one- 
dimensional space-fractional Gross-Pitaevskii equation can be defined as 

th^^ (x, t) = (-|^V^ + V{x)+g\tlj (x, tf\\ ^ (x, t) . (17) 
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Finally we can write time and space-fractional Gross-Pitaevskii equation in the one- 
dimensional form 



D]-'' (-^"^^ + V{x)+g\ij (x, tf\\ ^ (x, t) 



ifiip (x, t) = 

Here we focus solutions of the time-fractional Gross-Pitaevskii equation for different 
potentials. But we will consider other cases of these equations for different potential in 
the another study. 

5. Basic definitions for fractional calculus 

Before we discuss solution of the time-fractional Gross-Pitaevskii equation (16) we will 
give basic rules for calculating fractional differential equations in this section. Some 
definitions for fractional calculus are given below [22| H3| HH HHl HHl HTJ B8] . 

Definition 1 The Riemann-Liouville fractional integral operator of order a > of 
the function f{t) G C^ is defined as for /i > —1 

^oDrfit) = T;^j'dt'-^^, a>0, t>0 (19) 

i [a) Jto [t — t') 

where tgD^^'f (t) is the Riemann-Liouville fractional integral operator which is a direct 
extension of Cauchy's multiple integral for arbitrary complex a with Re (a) > 0. 
A fractional derivative is then established via a fractional integration and successive 
ordinary differential according to 

dt 

with Re (/3) > 0, and the natural number n satisfies the inequality n > Re (/3) > n — 1. 
Two special cases are the Riemann-Liouville o-Df for tg = and Weyl operator -oo-Df 
for to = — oo. The Riemann-Liouville operator for to = is given by 

1 (i*^ /■* f (f) 

oD^f (t) = — -^ / dt' ■' \\ „ , n > Reia) >n-l. (21) 

For any function / (x, t) the fractional Riemann-Liouville differintegration oDl~° = 
{d/dt)^ D^" is defined through the relation 



,,D'U it) = ^Dty it) (20) 



This differintegration operator of an arbitrary power for to = is given by 

which coincides with the heuristic generalization of the standard differentiation 
at^ [m — n)\ 



oA'a = — -t-^ . (25) 
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by introduction of Gamma function. An interesting consequence of Eq. (22) is the 
nonvanishing fractional differintegration of a constant a 

1 
T {a — u 
The Riemann-Liouville differentiation of exponential function leads to 

oA'e* = ^J^^^ Fi,! (1, l-u,t) (26) 

involving the confluent hypergeometric function Fi i. 

Definition 2 The Caputo fractional derivative operator [19] is given by 

I ft f (") (f) 

''^'l^'^'W^Jj'jhr^' n-l<a<« (27) 

Under natural condition on the function / (t), for a ^ n the Caputo derivative becomes 
a conventional ?7,-th derivative of the function / (t). 

The Riemann-Liouville fractional operator (22) and Caputo's fractional operators 
(27) have different form. Another difference between these operators is that the Caputo 
derivative of a constant is 0, whereas the Riemann-Liouville fractional derivative of a 
constant is given by Eq. (25). On the other hand, the main advantage of Caputo's 
approach is that the initial conditions for fractional differential equations with Caputo 
derivatives take some form as for integer-order differential equations. 

Definition 3 The Mittag-Leffler function [501 EI] -^a {z) is a complex function which 
depends on on two complex parameters a and /3. It may be defined by the following 
series representation when a > 0, valid in the whole complex plane 

E.A^)=]:^J^;^y «>0, ..C (28) 

In the case a and /3 are real and positive, the series converges for all values of the 
argument z, so the Mittag-Leffler function is an entire function. Some special cases 
of the Mittag-Leffler function are follow. The Mittag-Leffler function is the natural 
generalization of the exponential function. Being a special case of the Fox function 
introduced below, it is defined through the inverse Laplace transform 

Ea{-{t/Tf) = £-^l ^^1 (29) 

from which the series expansion 

can be deduced, where r is a constant. The asymptotic behavior of the Mittag-Leffler 
function interpolates between the initial stretched exponential form 

^<.(-('An~-p(-ff^) (31) 
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for t < T and the long-time inverse power-law behavior as 

B„ (_(,/,„. (i)°_L_ (32) 

for t^r, 0<a<l. Special cases of the Mittag-Leffler function are the exponential 
function 

E, i-t/r) = e-*/- (33) 

and the product of the exponential and the complementary error function 

Ei/2 (- it/rf') = e-'/^erfc ({t/rf') . (34) 

6. Solution of Gross-Pitaevskii equations of integer and fractional order 
with HPM 

The principles of HPM and its applicability for various kinds of differential equations 
ar given in|2a[25l[26l[27l|28l[29l[30l[33l[Ml[35l[Mll371|3^ Here after 

we present a review of the standard HPM and modified HPM suggested by Momani 
and Odibat pH [25]. We will employe HPM to the Gross-Pitaevskii equations of integer 
order and fractional order. 

6.1. HPM for Gross-Pitaevskii equation of integer order 

To explain the basic ideas of this method on nonlinear differential equation of integer 
order, we consider one-dimensional Gross-Pitaevskii equation in the form Eq. (1) and 
adopt the homotopy perturbation method to this equation. Gross-Pitaevskii equation 
(1) is a nonlinear partial differential equation which can be decomposed as linear and 
nonlinear part: 

L{^{x,t)) + N{^{x,t)) = f{x,t), xen, t>0 (35) 

where L{iIj) is a linear part and N{tp) is a nonlinear part and on the other hand / {x, t) 

is known analytical function of Eq. (1). For / {x,t) = 0, L{iIj) and N{iIj) can be defined 

follow as 

dib (x t) Ti^ 

L (V>) = ^^^^^^, N (^) = +^VV {x, t)+V (x) ^ (x, t)+g |^ (x, tf\ tfj {x, t) .(36) 

We note that Eq. (35) must be defined with the boundary condition 

B{i;{x), dip{x)/dn) =0, x eV (37) 

where F is the boundary of the domain Q. Homotopy perturbation method defines the 
homotopy (j){x,p) : Q x [0, 1] — )■ i? which satisfies 

H{<P,p) = il-p) [L(0) -L(^o)] +p[L(0) + iV(0) -fix,t)] = (38) 

or 

H (0, p) = L{<P)-L (^o) + VL (^o) + p [iV (0) - / (x, t)] = (39) 
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where p G [0, 1] is an embedding parameter, ipQ = ipi^x, 0) is an initial guess of exact 
solution ?/'(x,t), which is independent of time and satisfies the boundary conditions. 
Hence from Eqs. (38) and (39) we obtain 

i7(0,O) = L(0)-L(V^o(x)) = O (40) 

H{<P, 1) = L(V^o) + iV(0) - /(x, t) = . (41) 

Changing process p from zero to unity is just that of (p{x,p) from ipoi^x) to 0(x,t). 
In topology, this is known as homotopy [211 [251 ISSl |27], L{(j)) — L{iPq{x)) and 
L{ipo) + N{(f)) — f{x,t) are homotopic. This basic assumption is that the solution of 
Eqs. (38) or (39) can be expressed as a power series in p: 

(t){x,p) =(t)o{x) +p(t)i{x) +p^(t)2{x) + ... (42) 

The approximation solution of Eqs. (1) or (35) can be obtained in the limit p — )■ 1 as 

i'i.x.t) = (l){x,p)p__^^ = (1)0 (x) +p(j)i{x) +p^02(a;) + ••• • (43) 

As a consequence, the final result of standard Gross-Pitaevskii equation in the form 
Eq. (1) is given 

oo 

^P{x,t) = J2P^<pJi^) (44) 

j=o 

when p — 7- 1. Here we note that the solution function (p^ (x) includes the constants of 
the Eq. (1). On the other hand, we say that the convergence of the series (44) has been 
proved for many nonlinear ordinary and partial differential equations [26l ^^ . 

After a brief introduce of HPM, now we can apply HPM to the Gross-Pitaevskii 
equation of integer order. To solve Eq. (1), using Eq. (38) we construct the following 
homotopy: 

^-V^(p-V{x)(l)-g 



dcj) di)Q\ 
^^-^^[-d-t-^r 



0.(45) 



dt h \2m 

Substituting Eq. (42) into Eq. (45) and equating the coefficients of the terms with the 
identical powers of p, 

_ 900 dipo 



p 



p 



p 



dt 



dt 



0, 00 {x, 0) = ipo (x) = ip (x, 0) 



h 



1-R 2^^''^»-^'W^»-«Kk» 



902 i ( h 2, T./ N 



p 



3 . 



dt 



i_ Ji_ 
h V 2m 



V20o - V (x) 



i=0 k=0 

2 2-i 

^2-^EEI'^^ 

i=0 A;=0 



n-k-i 



n-k-i 



(46) 
(47) 

(48) 
(49) 



(50) 
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^' 4 ( ^v^<^,- - V (x) 0,_, - /e e" !</> 



dt h \ 2m 



'i\ \<Pk\ <Pj-k~i~l 
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= 0(51) 



i=0 fc=0 



(52) 



we get the iterative equation 



.i-V {x)(f)j-i- gJ2 J2 \<Pi\\<Pk\<Pj-k- 

i=0 k=0 



Adt, J = 1,2, 3,. ..(53) 



with initial value 

(pj {x, 0) = ^ (x) = ^ (x, 0) . (54) 

Hence all (pj values can be obtained using Eq. (53). The final solution of Gross-Pitaevskii 
equation (1) of integer order can be written in terms of the Eq. (53) in the limit p — ?> 1: 

rf / a 



*(-.*) = ^E/y„(^vv.-^w 



j=0 fc=0 



'k\(j)j-k-i dt .{55) 



HPM supposes that this solution satisfy Eq. (1) in the limit of the j9 — )■ 1. 



6.2. HPM for time-fractional Gross-Pitaevskii equation 

To solve the time-fractional Gross-Pitaevskii equation (16), we will apply HPM to 
Eq. (16). Hence we can construct the following homotopy 

? / h 



^^t<l>-j-[^'^'<l>-V{x)<j>-gr 



(l-p)(oA>- oA"^o)+P 
Suppose the solution of Eq. (56) to be as following form 

r {x,p) = 0? {x)+pri {x)+p'<p^ (x) + ... . (57) 

In the limit p — )■ 1 we can write solution of the time- fractional equation (16) 

oo 

r{x,t) = J2p^<P]{x) . (58) 

i=o 

Substituting Eq. (57) into Eq. (56) and equating the coefficients of the terms with 
identical powers of p, we obtain following equations 



0(56) 



n.0 . _ n"^a 



P" ■■ oDt 



dt 







p' : oA>? - { (^V^*^^ - \/ (x) 0^ - (? \cPl\ 0o) = 



1 l-i 



p": oD^>^- 1 ( Av20i-K(x) 01 -(?^^ 10,1 10,10 



1-fc-j 



(59) 
(60) 

0(61) 
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p' : oA>3 - 'j- (^V'<^2 -V{x)(l>,-gJ2J2\<t>^\ \<f>k\ 01-fc-.) = (62) 



(63) 



-'' ^'^ EE 

j=0 fc=0 



p': oDJ^0J-^|--VVi-^(^)</'.-i-^E E \H\^k\h~k~i\ =0(64) 



(65) 
Finally we can write iterative relation as 



j\^-^y (^VVi-^(^)0.-1-^E'E \<P^\\<Pk\<P,-k-l\dT{m) 



' ^r (a) JO \^z7u — ^0 J 

with initial value 

0J(a;,O) = ^(x)=^(a;,O), j = 1,2,3,.... (67) 

The final solution of time-fractional Gross-Pitaevskii equation (16) can be given in terms 
of the Eq. (66) in the limit p — )■ 1: 

CO 1.1 / t j-ij-i-i 



'ii i«; j=o -^^ \^'" i=o fc=o 

HPM supposes that this solution satisfy the Eq. (16) in the limit of the p — ?■ 1. 



7. Numerical implementations for different potentials 

Now we will investigate GP end time-fractional GP equations to discuss the ground 
state time-dependent dynamic of the Bose-Einstein condensation of weakly interacting 
system which involve attracting and repulsive interactions for different potential. In 
analysis we will use HPM to solve GP equations. 

Example 1 Firstly we consider GP equation of integer and fractional order for 
V {x) = and repulsive interaction g = 1 > 0. Hence we will discuss GPE of integer 
order in 1. Case- 1 and GPE of fractional order in 1. Case- 11. 

l.Case-I For V {x) = and repulsive interaction g = 1 > 0, the iterative equation (53) 
is given as 

i j-t f h :i-lj-^-l 

/o 



i ft % j-^j-i-^ \ 



i=0 fc=0 
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For g = 1 > 0, the HPM solutions of (pj are given by 

00 = tanhx (70) 

—it 

tanhx (71) 



tanhx (72) 

tanhx (73) 



r(2) 


t' 


r(3) 


tt' 



r(4) 



(74) 



(itV 
0n = (-!)" Trf^tanhx (75) 

1 (n + 1) 

where we set h = m = 1 for simphcity. If we put these terms in the series expansion 

oo 

^ {x, t) = 00 (x) + P01 (x) + p^02 (x) + ... = ^ pVi (a;) (76) 

i=o 

we obtain 

°^ „ (it)" 

■0„>o (x,t) = tanhx + T^ (—1)" — 7 tanhx . (77) 

^1 r (n + 1) 

This solution can be represented in terms of Mittag-Leffler function for a = 1 as 

■0g>o (x, t) = -El — ( — ) tanhx = tanhxe^*'^ (78) 

where i = l/r. This unnormalized solution is of the same analytical form with Eq. (11) 
which corresponds dark soliton solution of GP equation. The numerical demonstration 
of the Eq. (78) fot t = and the probability density of ipg^o {x, t) are given is given in 
Figs, (la) and (lb), respectively. As it can be seen from Fig. (lb) when external potential 
V (x) = 0, the probability density for attracting interactions has dark soliton form. 
These results confirm that obtained solution using HPM is the same with analytical 
solution of GP equation of integer order for V (x) = and (? > 0. 

l.Case-II Now we consider time-fractional GP equation for V (x) = and repulsive 
interaction i.e. g > 0. The iterative equation (66) of the time-fractional GP equation 
for V (x) = and (yf = 1 > is given by 

In this case, the HPM solutions of 0" are given by 



] 



tanhx (80) 
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Figure 1. In the case V {x) — and g > a) Numerical representation of the 
'0g>o (s^j 0)i b) Numerical representation of the probability density of '0g>o {x, t). 



if 



r(a + l) 



1.2a 



r(2a + i; 



tanhx 



tanhx 



it 



3o 



r (3a + 1) 



tanhx 



(81) 
(82) 
(83) 



lA) 



-1)"" -zr-, — tanhx 



5) 



r [na + 1 
where we set ^ = m = 1 for simplicity. If we put these terms in the series expansion 

(86) 



r (x,t) = 0° (x) +P0? {x) +p'<f)^ (x) + ... = ^p^0j (x) 

we obtain final result 

(zr)" 



V',%(x,t)=tanhx + 5](-l) 



n=l 



r (na + 11 



tanh X . 



^7) 



This solution can be written in terms of Mittag-Leffler function for arbitrary a value as 

CO / /j. / \o:\n 

^° (x,t) = tanh x + y \^ , ' \ tanh a; = E„ (- (t/r)°) tanh x (88) 

;^i r [na + 1) 

where we set i = 1/r". According to Eqs. (31) and (32) the asymptotic behavior of the 
Mittag-Leffler function, the fractional solution (88) is given by stretched exponential 
form 

(t/rT 



^g>o i^^ ^) = Ea (- (t/^T) tanhx ~ exp I p / ^ ^. 
or for long-time regime is given by inverse power-law as 



tanh a; 



/t\°' 1 
'^9>o i^^t) = Ea (- (t/r)°)tanhx ~ ( -j — — — — tanhx . 



(89) 



(90) 
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The nice analytical results are solutions of the time-fractional GP equation for V [x) = 
and g = 1 > 0. It can be clearly seen from Eqs. (89) and (90) that the time-dependent 
solution of time-fractional GP equation (16) is different from standard GP equation (1). 
Indeed, solutions (89) and (90) of time-fractional GP equation indicate that the ground 
state dynamics of the Bose-Einstein condensation evolute with time in complex space 
as obey to stretched exponential for short time regime and power law for long-time 
regime in the case external potential zero and interactions between bosonic particles are 
repulsive. Whereas for the same case i.e. V (x) = and g = 1 > the ground state 
dynamics of the condensation exponentially evolutes with time. However, the spatial 
part of both solution in Eq. (88) and in Eq. (78) are equal. Therefore, it is excepted 
that the time-fractional dynamics of condensation also produces dark soliton behavior 
similar to Eq. (78) as it can be seen in Fig. (lb). On the other hand, here we must 
remark that the fractional parameter a is a measure of fractality in between particles in 
condensation progress. Hence it determine time evolution of the ground state dynamics 
of condensation. 

In order to demonstrate the effect of the fractional parameter on the ground state 
dynamics or another say to investigate the ground state dynamics of the fractional 
condensation we plot the real and imaginer part of the ipgyo (^) ^) ^^^ several a values. 
As it can be seen from Figure 2 that the time evolution of the real and imaginer part 
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Figure 2. In the case V (x) — and 5 > a) Real part of the V'g>o i^^ ^) ^^r 
a = 1, 0.9, 0.8 values, b) Imaginer part of the V'o>o i^^ ^) ^o^' ^^ ~ Ij 0-9j 0-^ values. 



of the wave function ^g>o (^' ^) clearly depend on the fractional parameter a. For small 
values of the time, wave function solution depend on time nearly coincident for different 
a values, however fractional parameter a values substantially affect the solution of 
'^g>o (^' ^) when time is increased. Indeed for a < 1 values this effect can be clearly seen 
in Figure 2 (a) and (b) at large time values. 

On the other hand to see effect of the fractional parameter on the ground state 
dynamics of condensation of bosonic particles with fractality, we give the surface plot 
in Figure 3. This figure clearly shows that the real part of the wave function '?/'°>o i^i ^) 
evolutes to make oscillation with time for arbitrary a values, however, for in the limit 
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Re[l/'g>o 



Figure 3. Time and fractional parameter a dependence of the real part of '0g>o i^^ t) 
for V {x) ^Q and g > 0. 



a — )■ 0, at the same time, or in the hmit t — t- the wave function V'g>o (^)^) goes to 
stationary state which means that the wave function in these hmit does not change. 
Hence in these hmit values, stationary solution of the time-fractional GP equation is 
independent of time. For this reason, the soliton solution of time fractional GP equation 
(16) in these limit is stable similar to soliton solution of GP equation (1). The same 
behavior appears in imaginer part of the the wave function V'g>o (^' ^) • 

Example 2 Secondly we consider GP equation of integer and fractional order for 
y (x) = and attractive interaction g = —1 <Q. Hence we will discuss GPE of integer 
order in 2. Case-I and GPE of fractional order in 2. Case-II. 

2.Case-I For V (x) = and attractive interaction g = —1 < 0, the iterative equation 
(53) is given as 

"^^ = W I ^VVi -/e'E \<P^\ \<Pk\<P,-k-i] dt . (91) 

n JO y/m .^g ^^^ y 

In this case, the HPM solutions of (pj are given by 

00 = sec hx (92) 

it 
(hi = -— sec hx (93) 

sec hx (94) 



3r(3) 
it^ 



,^ , ,, sechx (95) 

4r (4) ^ ^ 



(96) 
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-1) 



sec hx 



where we set % 



m 



n + l)V{n + l) 
1 for simplicity. Using the series expansion 



i) (X, t) = 00 (x) + pcpi (x) + p^02 (x) + ... =J2 P'^Pj 

j=0 

the final results can be written as 

{^tr 



^g<0 {X, t) 



sec 



hx+Y.{-ir 



sechx 
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(97) 

(98) 
(99) 



^g^0{X,t)=E, 



sec hx = sec hxe *' ^ 



(100) 



in + l)Tin + l] 
This solution can be defined Mittag-Leffler function for a = 1 as 

where i = 1/r. This unnormalized solution is of the same analytical form with Eq. (13) 
which corresponds bright soliton solution of GP equation. The numerical demonstration 
of the Eq. (100) for t = and the probability density of ijjg^Q (x, t) are given is given in 
Figs. (4a) and (4b), respectively. As it can be seen from Fig. (4b) when external potential 
V (x) = 0, the probability density for attractive interactions has bright soliton form. 
These results confirm that obtained solution using HPM is the same with analytical 
solution of GP equation of integer order for V {x) = Q and (7 < 0. 





Figure 4. In the case V {x) — and g < a) Numerical representation of the 
'0g<o [x, 0), b) Numerical representation of the probability density of '0g<o (2;, t). 



2. Case-II Now we consider time-fractional GP equation for V (x) = and g = —1 < 0, 
The iterative equation (66) of the time-fractional GP equation for V (x) = 



Jo 



T 



\a-l 



i=0 fe=0 



-' hV (a) Jo ^ ' \ 2m 
For attractive interaction, the HPM solutions of 0° are given by 
00 = sec hx 
it 



'r(a + 2) 



sec hx 



10^1 0tfc-i)rf^ -(101) 

(102) 
(103) 
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^2a 



sechx (104) 



r (2a + 2) 

sechx (105) 



it^" 



r(3a 



(106) 



where we set /i = m = 1 for simplicity. All solutions can be added in series expansion 
as 

oo 

r {x,t) = 0^ (x) +pr, (x) +p'<p^ (x) + ... = ^p^0j (x) . (108) 

This produces the final results for g < 

r,<o {x^ t) = sec hx+J2 (-1)" rf \ o\ ^^" ^^ " (^O^) 

^1 r [na + 2) 

This solution can be written in terms of Mittag-Leffier function as 

oo / /i / \a\'n- 

C<o (x, t) = sec hx+J2 7 ^ A(./ ' , -,, sec /ix = E^ (- (t/r)") sec hx{nO) 

^ ^1 (na + l)r [na + 1) 

where we set i = 1/r". According to Eqs. (31) and (32) the asymptotic behavior of the 
Mittag-Leffler function, the fractional solution (110) is given by stretched exponential 
form 

r(« + i), 

or for long-time regime is given by inverse power-law as 

<<o {x, t) = E^ (- (t/r)") sec hx ~ (^^) t (I - a) ^'''' ^"^ ' ^^^^^ 

The nice analytical results are solutions of the time-fractional GP equation for V (x) = 
and g = —1 < 0. It can be clearly seen from Eqs. (Ill) and (112) that the time- 
dependent solution of time-fractional GP equation (16) is different from standard GP 
equation (1). Indeed, solutions (HI) and (112) of time-fractional GP equation indicate 
that the ground state dynamics of the Bose-Einstein condensation evolute with time in 
complex space as obey to stretched exponential for short time regime and power law for 
long-time regime in the case external potential zero and interactions between bosonic 
particles are attractive. Whereas for the same case i.e. V (x) = and g = —1 < the 
ground state dynamics of the condensation exponentially evolutes with time. However, 
the spatial part of both solution in Eq. (110) and in Eq. (100) are equal. Therefore, it is 
excepted that the time-fractional dynamics of condensation also produces bright soliton 
behavior similar to Eq. (100) as it can be seen in Fig. (4b). On the other hand, here 



^^^Ja:,t) ^E„(-(t/r)°)sec/ia;~exp {- ^)[ ' ^, \sechx (111) 
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we must remark that the fractional parameter a is a measure of fractahty in between 
particles in condensation progress. Hence it determine time evolution of the ground 
state dynamics of condensation. 
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Figure 5. In the case V {x) — and g < a) Real part of the ■ipg^Q{x,t) for 
a = 1, 0.9, 0.8 values, b) Imaginer part of the V'g<o i^^ ^) fo^' ct = 1; 0-9j 0-8 values. 



In order to demonstrate the effect of the fractional parameter on the ground state 
dynamics or another say to investigate ground state dynamics of fractional condensation 
we plot the real and imaginer part of the V'g<o (^)^) fo^ several a values. As it can be 
seen from Figure 5 that the time evolution of the real and imaginer part of the wave 
function V'g<o i^^t) clearly depend on fractional parameter a. For small values of the 
time, wave function solution depend on time nearly coincident for different a values, 
however fractional parameter a values substantially affect the solution of V^g<o l-^' ^) when 
time is increased. Indeed for a < 1 values this effect can be clearly seen in Figure 5 (a) 
and (b) at large time values. Furthermore, to see effect of the fractional parameter on 




Re[^/Co] 



Figure 6. Time and fractional parameter a dependence of the real part of ^g<o i^: t) 
for y (cc) = and g <Q. 
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the ground state dynamics of condensation of bosonic particles, we give the surface plot 
in Figure 6. This figure clearly shows that the real part of the wave function ipg^Q {x, t) 
evolutes to make oscillation with time for arbitrary a values, however, for in the limit 
a — )■ 0, at the same time, or in the limit t ^ the wave function ipg^Q {x,t) goes to 
stationary state which means that the wave function in these limit does not change. 
Hence in these limit values, stationary solution of the time-fractional GP equation is 
independent of time. For this reason, the soliton solution of time fractional GP equation 
(16) in these limit is stable similar to soliton solution of GP equation (1). The same 
behavior appears in imaginer part of the the wave function ipg^Q {x,t). 

Example 3 Thirdly, in this example, by using HPM we will obtain solutions of the 
GP equation of integer and fractional order for optical lattice potential V {x) = sin^ x 
and repulsive interaction i.e. g = 1 > 0. Hence we will consider GPE of integer order 
in 3. Case-I and GPE of fractional order in 3. Case-II. We note that we have found the 
analytical solution of GP equation as il!g=i (a;, t) = e~^**/^ cos a; in the case 5^ = 1 > for 
optical lattice potential V {x) = sin^x which satisfy Eq. (1). 

3. Case-I For optical lattice potential V (x) = sin^ x and repulsive interaction g = 1 > 
0, the iterative equation (53) is given as 




i\ \(Pk\ (Pj-k-1 



For this potential, the HPM solutions of (pj are given by 

00 = cosx (114) 

Sit 

01 = -^pT^cosa; (115) 

cos a; (116) 



4r(3) 
27it^ 



8r(4) 



cosa; (117) 



11^ 



0„= -1 " ^ ' -cosa: 119 

where we set /i = m = 1 for simplicity. If these 0j terms can be put in the series 

00 
^ (a;,t) = 00 (a;) +P01 (x) +^^02 (x) + ... = ^p^i {^) (120) 
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(3it)" 



■ cos X . 



we obtain final results for g = 1 > 

oo 

■000 (2^, ^) = cosx + y^ (— 1)' ^, 
The solution (121) can be written in terms of Mittag-LefHer function for a = 1 as 

(tMr 



:i2ii 



^g>0 {X, t) =COSX+J2 



n=l 



r(n + i] 



■ COS X = Ei{— {t/r)) cos x 



(122) 



where we set 3i/2 = 1/r. Eq. (122) can be given in exponential form as 

^g>o (a;, t) = El (- (t/r)) cos x = e"*/^ cos a; . (123) 

Our analytical result in Eq. (123) is the same with analytical solution of GP equation 
for optical lattice potential V (x) = sin^ x and attractive interactions g > 0. The 
numerical demonstrations of the Eq. (123) for t = and the probability density of 
ipg^o {x,t) are given is given in Figs. (7a) and (7b), respectively. As it can be seen from 
Fig. (7a) that the solution ^g>o {x, t) for t = has a Gaussian form, whereas probability 
density of V'g>o i^, t) has a double well shape for optical lattice potential V (x) = sin^ x 
and attractive interactions g = 1 > 0. As a result we note that these analytical and 
numerical results confirm that obtained solution using HPM is the same with analytical 
solution of GP equation of integer order for V (x) = sin^ x and g = 1 > 0. 
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Figure 7. V (x) — sin^ x and repulsive interactions g = 1 > 0, a) Numerical 
representation of the ?/'g>o (x, 0) , b) Numerical representation of the probability density 
oiiljg>o{x,t). 



3. Case-II Now we consider time-fractional GP equation for V (x) = sin^ x and 
repulsive interaction g = 1 > 0. The iterative equation (66) of the time-fractional 
GP equation for this potential is given by 



it 



T 



\a-l 



h 



^^ hT (a) Jo '^ ' ' \1m 
For (yf = 1 > 0, the HPM solutions of 0^ are given by 

0Q = COS X 



V'<PJ_i - sin^x^J^, - ^E E ICI 10^1 <PU^i dr (124) 



i=0 k=0 



(125) 
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cosx (126) 

-cos a; (127) 



2r(a + l) 



4r (2a + 1 
27it3° 



cos a; (128) 



Sr (3a + 1) 

(129) 



2"r (na + 1' 



■cos a; (130) 



where we set h = m = 1 for simplicity. These terms can be put in a seises to obtain 

^"(a;,t) 

oo 

r {x,t) = 0^ (x) +pr, (x) +p'<p^ (x) + ... = ^pvj (x) (131) 

i=o 
We obtain final results foTg = l>0 

^" (x, t) = cosx + y (-1)" — ^-^ COS a; . (132) 

'^5>ov , J ^1 2«r(na + l) ^ ' 

This solution can be written in terms of Mittag-Leffler function as 

ipa>() {x,t) = COS a; + T^ — cos a; = Ea (— (t/T)'^) cos a; (133) 

^1 r [na + 1) 

where we set 3i/2 = 1/r". According to Eqs. (31) and (32) the asymptotic behavior 
of the Mittag-Leffler function, the fractional solution (133) is given by stretched 
exponential form 

r(a + l), 

or for long-time regime is given by inverse power-law as 

/t \ " 1 
■^ g>o{x, t) = Ea{-{t/Tf) cosx ^ i-j pQ_^x cosx ■ (135) 

Eqs. (134) and (135) are solutions of the time- fractional GP equation for V {x) = sin^ x 
and (^ = 1 > 0. It can be clearly seen from Eqs. (134) and (135) that the time-dependent 
solution of time-fractional GP equation (16) is different from standard GP equation 
(1). Indeed, solutions (134) and (135) of time-fractional GP equation indicate that for 
a < 1 the ground state dynamics of the Bose-Einstein condensation evolute with time in 
complex space as obey to stretched exponential for short time regime and power law for 
long-time regime in the case external potential V (x) = sin^ x and interactions between 
bosonic particles are repulsive. Whereas for the a = 1, the ground state dynamics 
of the condensation exponentially evolutes with time. However, the spatial part of 



^^^g (x,t) = E^ (- (t/r)") cos a; ~ exp ( - ^;; ; ^^ ) cos a; (134) 



Solution of the time-fractional Gross- Pitaevskii equation 



22 



both solution in Eq. (133) and in Eq. (123) are equal. Therefore, it is excepted that 
the time-fractional dynamics of condensation also produces double well shape behavior 
similar to Eq. (135) as it can be seen in Fig. (7b). On the other hand, here we must 
remark that the fractional parameter a is a measure of fractality in between particles in 
condensation progress. Hence it determine time evolution of the ground state dynamics 
of condensation. In order to demonstrate effect of the fractional parameter on the ground 





(a) 



(b) 



Figure 8. In the case V (x) — sin x and g > a.) Real part of tlie '0g>o (^i f^i' 
a = 1, 0.9, 0.8 values, b) Imaginer part of the V'o>o (^' ^) ^^^ a — 1, 0.9, 0.8 values. 



state dynamics or another say to investigate ground state dynamics of condensation with 
fractality we plot the real and imaginer part of the V'o>o (^' ^) ^"^^ several a values. As it 
can be seen from Figure 8 that the time evolution of the real and imaginer part of the 
wave function ^g>o {x,t) clearly depend on fractional parameter a. For small values of 
the time, wave function solution depend on time nearly coincident for different a values, 
however fractional parameter a values substantially affect the solution of ipg-^Q {x, t) when 
time is increased. Indeed for a < 1 values this effect can be clearly seen in Figure 8 




ReW" 



Figure 9. Time and fractional parameter a dependence of the real part of '0g>o i^^ ^) 
for V {x) = sin^ x and g > 0- 



Solution of the time-fractional Gross- Pitaevskii equation 23 

(a) and (b) at large time values. Furthermore, to see effect of the fractional parameter 
on the ground state dynamics of condensation of bosonic particles with fractality, we 
give the surface plot in Figure 9. This figure clearly shows that the real part of the 
wave function ipg^Q {x, t) evolutes to make oscillation with time for arbitrary a values, 
however, for in the limit a — )■ 0, at the same time, or in the limit t — )■ the wave function 
ipgyo {x, t) goes to stationary state which means that the wave function in these limit 
does not change. Hence in these limit values, stationary solution of the time-fractional 
GP equation is independent of time. 

Example 4 Finally, we will obtain solutions of the GP equation of integer and 
fractional order for optical lattice potential V (x) = — sin^ x and attractive interaction 
i.e. g = —1 < 0. Hence we will consider GPE of integer order in 4-Case-I and GPE of 
fractional order in ^.Case-II. We note that we have found the analytical solution of GP 
equation as ^g=-i {x, t) = e**/^ cos a; in the case g = —1 < for optical lattice potential 
V (x) = — sin^a; which satisfy Eq. (1). 

4.Case-I For optical lattice potential V (x) = — sin^x and attractive interaction 
g = —1 < 0, the iterative relation (53) is written as 

i ft f h j-'i^j-i-'^ \ 

(t>3 = T ^VVj-i + sin^ x(f)j_i -gY^ Y^ |0i| |0fe| <Pj-k-i dr (136) 

tiJo yzm .^Q ^^g J 

For this potential, the HPM solutions of (pj are given by 

00 = cosx (137) 

it 

01 = ————— cosx (138) 



2r(i) 
t^ 

'at (3) 
it^ 

'8r(4) 



cosx (139) 

cos a; (140) 



(141) 



' 2"r(n + ll 



cosx (142) 

where we set ^ = m = 1 for simplicity. If these terms are put in the series 



i){x,t) = 00 (x) +p0i(x) +p^02(x) + ... = Yp'^^A^) (143) 

j=0 

we obtain final results as 

'0o<o ix,t) = cosx + y^ (—1)"" ; -cosx . (144) 
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The solution (144) can be written in terms of Mittag-Leffler function for a = 1 as 

(t/r)) 



1pg<oix,t) = COSX+ ^ 



cosx = Ei{— (t/r)) cosx 



^1 r {na + 1) 

where we set i/2 = 1/r. Eq. (145) can be given in exponential form as 
fpg<o i^^t) = El (— {t/r)) cosx = e~*''^cosa: . 



(145) 



(146) 
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Figure 10. V {x) = — sin^ x and attractive interactions 5 = — 1 < 0, a) Numerical 
representation of the ipg<o (^^j 0) , b) Numerical representation of the probability density 



Our HPM result in Eq. (146) is the same with analytical solution of GP equation for 
optical lattice potential V {x) = — sin^x and attractive interactions g < 0. At t = 0, 
the numerical demonstrations of the ipg<o {x, 0) and the probability density oiipg<:o (x, 0) 
are given is given in Figs. (10a) and (10b), respectively. As it can be seen from Fig. (7a) 
that the solution 4'g<:o (x, 0) for t = has a Gaussian form, whereas probability density 
of ipg<:o (x, t) has a double well shape for optical lattice potential V^ (x) = — sin^ x and 
attractive interactions 5' = — 1<0. As a result we note that these analytical and 
numerical results confirm that obtained solution using HPM is the same with analytical 
solution of GP equation of integer order for V {x) = — sin^ x and g = —1 < 0. 



4.Case-II In last case, we will consider time-fractional GP equation for optical lattica 
potential V (x) = — sin^x and attractive interaction g = —1 < 0. The iterative equation 
(66) of the time-fractional GP equation for this potential is given by 

J^ it - rr-' [^^'<PU + sm^^0^-i - /e'e' \€\ m <PU^i ) dr{U7) 



6J 



hT (a 
For this potential, the HPM solutions of 0j are given by 

0Q = COSX 



i=0 A:=0 



2r(a + ll 



COSX 



i2o 



4r(2a + i; 



COSX 



(148) 
(149) 

(150) 
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it^"' 



cosx (151) 



^•^ 8r(3a + l) 

(152) 

C-(-ir^.r|^-- (153) 

where we set h = m = 1 for simplicity. If these terms are put in the series 

oo 
j=0 

we obtain final results as 

i^g<oix,t) = COS3: + ^(-l)" cosx . (155) 



n=l 



2«r {na + 1) 



V'g<0 (^5^) = -^a (~ (V''")") COSX ~ exp I — rp;-;^ ;— j-^ ) COS X (157) 



This solution can be written in terms of Mittag-Leffler function as 

oo / (f/'r\°'\'"' 

C<o i^^ t) = COS a: + V -— — — cos a; = E« (- (t/r)°) cos x (156) 

^1 r [na + 1) 

where we set i/2 = 1/r". According to Eqs. (31) and (32) the asymptotic behavior of the 
Mittag-Leffler function, the fractional solution (156) is given by stretched exponential 
form 

jt/rr 

r(a + l)^ 

or for long-time regime is given by inverse power-law as 

/t \ " 1 

-^gKoi^^t) = Ec.{-{t/ t)"") cos X r^ i-j pQ_^x cosx ■ (158) 

Eqs. (157) and (158) are solutions of the time-fractional GP equation for V (x) = — sin^ x 
and g = —1 < 0. It can be clearly seen from Eqs. (157) and (158) that the time- 
dependent solution of time-fractional GP equation (16) is different from solution of the 
standard GP equation (1). Indeed, solutions (157) and (158) of time-fractional GP 
equation indicate that the ground state dynamics of the Bose-Einstein condensation 
evolute with time in complex space as obey to stretched exponential for short time 
regime and power law for long-time regime in the case external potential V (x) = — sin^ x 
and interactions between bosonic particles are attractive. Whereas for the a = 1, the 
ground state dynamics of the condensation exponentially evolutes with time. However, 
the spatial part of both solution in Eq. (156) and in Eq. (146) are equal. Therefore, 
it is excepted that the time-fractional dynamics of condensation also produces double 
well shape behavior similar to Eq. (146) as it can be seen in Fig. (10b). On the other 
hand, here we must remark that the fractional parameter a is a measure of fractality 
in between particles in condensation progress. Hence it determine time evolution of the 
ground state dynamics of condensation. In order to demonstrate effect of the fractional 
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(a) 




t 
(b) 



Figure 11. In the case V (x) = — sin x and g > a) Real part of the ^g<o i^: t) for 
a = 1, 0.9, 0.8 values, b) Imaginer part of the V'g<o (^i t) for a — \, 0.9, 0.8 values. 




Re[<< 



g<0J 



Figure 12. Time and fractional parameter a dependence of the real part of V'q<o {x^ t) 



for V {x) 



sin X and 5 < 0. 



parameter on the ground state dynamics or another say to investigate ground state 
dynamics of condensation with fractahty we plot the real and imaginer part of the 
V'g<o (^5 ^) fo^ several a values. As it can be seen from Figure 11 that the time evolution 
of the real and imaginer part of the wave function '?/'"<o (^5 ^) clearly depend on fractional 
parameter a. For small values of the time, wave function solution depend on time nearly 
coincident for different a values, however fractional parameter a values substantially 
affect the solution of '?/'g<o {^i t) when time is increased. Indeed for a < 1 values this 
effect can be clearly seen in Figure 11 (a) and (b) at large time values. Furthermore, to 
see effect of the fractional parameter a on the ground state dynamics of condensation of 
bosonic particles, we give the surface plot in Figure 12. This figure clearly shows that 
the real part of the wave function ipg^o {x, t) evolutes to make oscillation with time for 
arbitrary a values, however, for in the limit a — >■ 0, at the same time, or in the limit 
t — )• the wave function i'g^^o {x, t) goes to stationary state which means that the wave 
function in these limit does not change. Hence in these limit values, stationary solution 
of the time-fractional GP equation is independent of time. 
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8. Conclusion 

In this study, we liave suggested fractional Gross-Pitaevskii equations to investigate the 
time-dependent ground state dynamics of the Bose-Einstein condensation of weakly 
interacting bosonic particle system which can includes non-Markovian processes or 
non-Gaussian distributions and long-range interactions. However we focused the time- 
fractional Gross-Pitaevskii equation and obtained solutions of the the Gross-Pitaevskii 
equation of the integer and fractional order for repulsive and attractive interactions in 
the case external trap potentials V{x) = and optical lattice potential V (x) = isin^a; 
using Homotopy Perturbation Method. Summarizing, in the Example 1 and 2, we have 
considered GP equation of integer and fractional order for the repulsive and attractive 
interactions in the case external trap potential V{x) = 0. On the other hand, in 
the Example 3 and 4, we have discussed GP equation of integer and fractional order 
for the repulsive and attractive interactions in the case the optical lattice potentials 
V (x) = ± sin^ X. We found several results: 1) we have shown that our HPM solutions of 
GP equation of integer order are the same the analytical solutions of its for V (x) = and 
optical lattice potentials V (x) = ± sin^ x. These confirms that HPM is reliable method 
to solve time-dependent non-linear partial differential equations as well GP equation, 
2) We show that the time dependent part of the solutions of the time-fractional GP 
equation have stretched exponential form for repulsive and attractive interactions in 
the case V {x) = and optical lattice potentials V (x) = ± sin^ x. However, for the 
same potentials and interactions, time dependent part of the solution of GP equation 
of integer order have exponential form. These means that time evolution of the ground 
state dynamics of condensation of bosonic particles including non-Markovian processes 
deviates exponential form, and evolutes with time as stretched exponentially. 

We conclude that the condensation of weakly interacting bosonic particles have non- 
Markovian processes which can be lead fractional dynamics in condensate phase. This 
dynamics can be represented by time-fractional Gross-Pitaevskii equation. Although 
time-fractional GP equation does not predict critical temperature Tc of the bosonic 
system, it can help to understand ground state dynamic of these system. To 
understand the contradictions of the experimental and theoretical predictions of critical 
temperature about condensation in real bosonic gas, the memory effects of the non- 
Markovian processes, non-Gaussian distribution of particles, the long-range inter- 
particle interaction effects or the fractal structure of the interacting in a bosonic system 
below the T^ can be taken into account. In the next study, to investigate of the 
condensation dynamics of weakly interacting bosonic systems we will consider the time 
and space-fractional GP equation for many different potentials, which have no analytical 
solutions. And also we will try to determine relation between condensation temperature 
of real bosonic system and a used in the fractional Gross-Pitaevskii equation. 
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